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ABSTRACT 

We report on a measurement of the temperature of the cosmic microwave background ra- 
diation field, Tcmb, at z = 0.88582 by imaging HC3N(3 «— 2) and (5 ■<— 4) absorption in the 
foreground galaxy of the gravitationally lens magnified radio source PKS 1830—211 using the 
Very Long Baseline Array and the phased Very Large Array. Low-resolution imaging of the 
data yields a value of T ro t = 5.6^7jg K, for the rotational temperature, T rot , which is con- 
sistent with the temperature of the cosmic microwave background at the absorber's redshift of 
2.73(1 + z) K. However, our high-resolution imaging reveals that the absorption peak position 
of the foreground gas is offset from the continuum peak position of the synchrotron radiation 
from PKS 1830—211 SW, which indicates that the absorbing cloud is covering only part of the 
emission from PKS 1830—211, rather than the entire core-jet region. This changes the line-to- 
continuum ratios, and we find T rot between 1.1 and 2.5 K, which is lower than the expected value. 
This shows that previous, T ro t, measurements could be biased due to unresolved structure. 

Subject headings: 



cosmic background radiation — cosmology: ob- 
servations — galaxies: ISM - quasars: absorption 
lines — quasars: individual (PKS 1830—211) 

1. INTRODUCTION 

PKS 1830-211 i s a well-studied high -redshift 



blazar (z = 2.507; iLidman et al~1ll999t) . which 



is gravitationally lensed by a foreground face-on 
spiral ga laxy revealed by Hubble Space Telescop 
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imaging (IWiklind k Combes ll99 8t lCourbin etaL 
200 alWinn et al. |2002h at z = 0.88582 (jWiklind k Com' 
19961) . The radio core of PKS 1830-211 is lensed 

into two components: the southwest (SW) and 

northeast (NE) images ([Pramesh Rao k Subrahmanvan 



1988 : Subrahmanvan et al. 1990f) with a separa 



tion of 0;'971 (|Jin et al. 1120031) or ~ 7.7 kpc at 
z = 0.89 (assuming a flat universe with the cos- 
mological parameters Hq = 70 km s _1 Mpc -1 , 
n m = 0.27 and D, A = 0.73). The two bright 
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core images are connected by a ring-like structure 
of low surface brightness that is often referre d 



1991), 



to as an Einstein ring ([Jauncev et al 
whose center is near the nucleus of the lensin. 

z 



0.89 galaxy (IKo chanek k Naravan I 11 992 
Nair. Narasimha. k Rao Ill993l : see Courbin et al. 
2002 for a historic summary). The SW and NE 
lensed images of PKS 1830—211 lie on both sides 
of the bulge of the lensing galaxy, with projected 
distances of - 3.0 kpc (0"40) a nd ~ 4.4 kpc (0"59 



fcym the nucleus, re spectively (jWinn et al. 11200 



nten et al. 



20081) . 

Radiation from the SW lensed component 
of PKS 1830-211 is strongly affected by the 
intervening gas in the foreground z = 0.89 
galaxy that causes reddening and strong ab- 
sorption features detected fo r a large number of 
different molecular species dWiklind k Combes I 
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Although the disk of the z = 0.89 galaxy 
is s een nearly face-on (incl ination i < 20°; 
e.g. JWikhnd k Combes1ll998l) . the velocity width 
of the absorption toward the SW image reaches 



1 



~ 100 km s- 1 FWZP (|Muller fc Guelin 1 120081) . 
which is significantly larger than that of most gi- 
ant molecular clouds (GMC) of the Milky Way 
and better resembles Galact i c Center clouds 
dMeriten. Carilli. fc Reidl 119991 iMuller fc Guelin 
2008). Toward the NE lensed image, weaker 
and narrower (~ 15 km s^ 1 FWHM) absorp- 
tion components have been detected at a velocity 
shifted by —147 km s _1 rel ative to the absorp- 



tion t oward the SW imag e (jWiklind fc Combes I 
19981: IMuller et al~l 12009) . The characteristics 



of the absorption in front of the NE image are 
similar to those of the Galactic diffuse clouds 



(IMuller fc Guelin 1 120081 ) 



High-redshift molecular absorption provides an 
important means to understand the physical con- 
ditions in the molecular interstellar medium in 
a galaxy at a cosmologically significant redshift. 
Of particular interest is an accurate measure- 
ment of the temperature of the cosmic microwave 
background (CMB) radiation (Tcmb), one of the 
fundamental cosmological pa rameters, and its 
varia tion with redshift (e.g., iNoterdaeme et al^ 
20111 ). The standard cosmological model pre- 



dicts, for an adiabatic expansion of the Universe, 
Tcmb(z) = TcMB-ofl+z), where Tqmb. o = 2.73 K 
(jMather et al. I Il994 Il999t iFixsen et al~l Il996f) . 
This important relation predicted by the theory 
needs to be confirmed by direct observations. This 
can be achieved by measuring the excitation tem- 
perature (T ex ) of interstellar atomic or molecular 
species at high redshift, in particular the rota- 
tional excitation temperature (T rot ) of molecular 
transitions. In diffuse interstellar gas, where col- 
lisional excitation is negligible, rotational tran- 
sitions can be in radiative equilibrium with the 
CMB (i.e., T rot = T CM b)- The rotational tem- 
perature T rot of the gas can be calculated by ob- 
serving two transitions of a given molecule, and 
therefore Tomb at the redshift of the gas can be 
obtained. For PKS 1830- 211, this was done by 
Wiklind fc Combes (Il996l) . who determined an 
upper limit of Tcmb < 6.0 K. They found that 
only 36% of the background continuum was cov- 
ered by the foreground absorbing gas, implying 
observations with higher angular resolution were 
needed to a chieve a reliable value for T ot . A few 
years later, ICombes fc Wiklind ] (fl999l) observed 
a larger sample of molecular transition lines, re- 
sulting in several measurements of Tcmb below 



the pr edicted value. Subsequently, Henkel et al. 
(2009) performed similar observations with the 
100-m Effelsberg radio telescope, deriving the 
covering factor and optical depth from advanced 
models fitted to the data. Their analysis was 
largely consistent with a Tcmb value of 5.14 K, al- 
though individual models co uld show large offsets. 
Finally, IMuller et al.l (1201 lh performed this anal- 
ysis with the Australia Telescope Compact Array, 
providing a higher resolution. Their findings were 
mostly consistent with Tcmb = 5.14 K, although 
some transitions showed significantly lower values 
than the theoretical prediction. Even higher reso- 
lution observations are needed to derive a precise 
measurement. 

In this paper, we report on a measurement of 
the temperature of the microwave background ra- 
diation field at z = 0.88582 by VLBI imaging of 
two microwave transitions of HC3N. The obser- 
vations and data reduction are described in Sec- 
tion [2] Our measurement of T rot is presented 
in Section [3j using both low- and high-resolution 
imaging (Sections 13.11 and 13. 2\ respectively), fol- 
lowed by a discussion of possible biases in Sec- 
tion 21 Finally, our conclusions are summarized in 
Section [3J 

2. OBSERVATIONS AND DATA ANAL- 
YSIS 

2.1. Observations 

We observed the redshifted HC 3 N(3 4- 2) and 
HC 3 N(5 <- 4) absorption lines at 14.5 GHz and 
24.1 GHz, respectively, toward PKS 1830-211 us- 
ing the National Radio Astronomy Observatory't0 
Very Long Baseline Array (VLBA) combined with 
the phased Very Large Array (VLA) as a VLBA 
antenna on 1999 January 29 (program BC087B). 
We employed two frequency bands of 8 MHz band- 
width each, in both right- and left-circular po- 
larizations. The center frequencies of the two 
bands were set to cover both velocities of the 
NE and SW components (with 14 7 km s _1 dif- 



ference; |Wiklind_jL£conbesJ [1998): respectively, 
14.46272 GHz and 14.47419 GHz for the J = 3 <- 
2 transition and 24.10444 GHz and 24.12356 GHz 
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for the J = 5 «— 4 transition. 

We used seven observing blocks of 1 hour dura- 
tion each. Subsequent observing blocks were un- 
dertaken at alternating frequencies of 14.5 GHz 
and 24.1 GHz (in total, four observing blocks at 
14.5 GHz and three at 24.1 GHz). Each 1 hour 
block consisted of five repetitions of 12-minute 
groups, each consisting of a 3-minute scan to- 
ward the quasar J1733— 1304 as passband calibra- 
tor and a 9-minute scan toward the target source 
PKS 1830-211. 

The data correlation was performed at the 
VLBA correlation facility in Socorro, NM. The 
data from each antenna pair were cross-correlated 
with an integration time of 2 seconds. Each 8 MHz 
band was split into 128 spectral channels, yield- 
ing a channel separation of 62.5 kHz. The corre- 
sponding velocity resolutions (with the relativistic 
Doppler effect) are 0.887 km s" 1 and 0.532 km s _1 
for the J = 3 <— 2 and J = 5 4 transitions, re- 
spectively, for rest fre quencies of 27294.289 MH z 
and 45490.3138 MHz (jMiiller et al. Il200ll [20051) . 
The data correlation process was repeated, adopt- 
ing two different positions for the correlation 
center of PKS 1830—211, i.e., the positions of 
the SW and NE images. In J2000 coordinates, 
the pointing position for the observation was 
(18 h 33 m 39?9,-21°03'40"5), and the adopted cor- 
relation positions were (18 h 33 m 39?8816, -21°03'40"6080) 
for the SW position and (18 h 33 m 39?9316, -21°03'39"8197) 
for the NE position. 

2.2. Calibration 

The correlated data were calibrated using the 
NRAO Astronomical Image Processing System 
(AIPS; lGreisenll2003h . We did not detect ab- 
sorption toward the NE component. Therefore, 
we analyzed only the data set with the correlation 
center at the SW component position and only the 
8MHz bands that covered the absorption spectrum 
in each of the 14.5 GHz and 24.1 GHz bands. 

We first corrected the interferometer delays and 
phases for the effects of diurnal feed rotation (par- 
allactic angle). We also removed the off-source 
VLA station data by flagging the first minute of 
each scan. The amplitudes of the interferometric 
data were corrected for the small biases (~ a few 
percent) of the data sampling thresholds among 
different stations, and then converted from corre- 



lation coefficients to flux density units using the 
system temperature and antenna gain curve infor- 
mation. 

We removed instrumental delay and phase off- 
sets by fitting fringe patterns to the data of the 
calibrator J1733— 1304 and corrected the data for 
these offsets. Bandpass amplitude and phase cor- 
rections were also determined from J1733— 1304 
and applied to PKS 1830-211. We then per- 
formed phase calibration by fitting fringe patterns 
to the data of PKS 1830—211 on time scales of 
20 seconds, using a model based on the image of 
PKS 1830-211 from a trial fringe fit. 

A continuum spectral baseline for PKS 1830—211 
was fitted based on channels 10 to 25 and 105 to 
120 and subtracted from the visibility data, yield- 
ing a continuum-only data set and a baseline- 
subtracted "line" data set. The phase of the 
continuum data set was further calibrated us- 
ing an iterative self-calibration process. During 
the self-calibration process, we simultaneously 
CLEANed the SW and NE continuum compo- 
nents. The NE component had a positional off- 
set of (Aacos<5, AS) = (0"6425, 0"7276) relative 
to the SW component and was still in the field 
of view of the data correlated at the SW posi- 
tion. The interferometer gain corrections from 
self-calibration were then applied to the baseline- 
subtracted "line" data set of PKS 1830-211. 

Before the spectrum fitting, we converted the 
frequency v of the observed spectra into the veloc- 
ity units v by using the relativistic Doppler effect, 
v = vq — w/c)/(l + v/c), where c is the speed 
of light and v$ the rest frequency of the transition. 
We subtracted the expected velocity of 168198.25 
km s _1 for the redshift of z = 0.88582, which de- 
fines the zero velocity of the obtained spectra. 



2.3. Rotational Temperature 

From the measured absorption amplitude AS(v) 
(defined to be negative for absorption) and contin- 
uum levels S c of the two transitions, HCaN(3 <— 2) 
and HCsN(5 <— 4), the optical depth t(v) of each 
transition is determined by 



t(v) = - In 1 + 



AS(v) 

fcS c 



(1) 



where f c is the covering factor of the background 
continuum source by the absorbing cloud ( f c < 1). 
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If we assume the two transitions are caused by the 
same absorbing cloud, we can use these values in 
the following equation to determine the rotational 
temperature T rot of the absorbing gas: 



-AEa 



1 



ln f33 93 r M (v ) (l-c-^32/^ cx) 
\35 95 r 32 (v ) (l-e-hv^/kT^ 



(2) 



where A.E42 is the energy difference between the 
,7 = 4 and J — 2 levels, k the Boltzmann's 



const ant (AE i2 /k = 3.06 K; iMuller et al. 112001 
2005), g.j = 2 J + 1 the statistical weight of the 
rotational level J, and v u i and t u i(vq) are the rest 
frequency and optical depth (at the line center ve- 
locity wo), respectively, of the rotational absorp- 
tion J = u <— I. We assume equal excitation tem- 
perature T ox for the two transitions J = 3 <— 2 and 
J = 5 <— 4 (i.e., local thermodynamical equilib- 



rium) and also assume T r , 



T ox since collisional 



excitation is negligible for the diffuse interstellar 
gas. We solved Equation [2] for T rot numerically 
through iteration. The equations used to derive 
Equation [2] are documented in the Appendix. 

3. RESULTS 

3.1. Low Resolution Analysis 

We first imaged the continuum and line absorp- 
tion with low spatial resolution in order to simu- 
late the results from previous studies. We made 
image cubes at both 14.5 GHz and 24.1 GHz with 
a circular restoring beam of 26 mas (FWHM); at 
this resolution, the continuum emission is essen- 
tially unresolved. Figure [1] displays these low res- 
olution images; the absorption distributions were 
obtained by plotting absorption depth, \AS\, at 
each pixel of the image. We performed Gaus- 
sian fits to the observed spectra of the baseline- 
subtracted "line" image cubes and measured the 
absorption amplitudes of the HCaN(3 <— 2) and 
HC3N(5 <— 4) transitions at each pixel of the im- 
age. 

We first solved for the three parameters, the 
amplitude, AS(v), the line width (FWHM), a v , 
and the central line velocity, vq, of the ab- 
sorption at the peak pixel of both continuum 
and absorption at the center of the map. We 
measured ^0,32 = —12.88 ± 0.32 km s -1 and 
(T V ,32 = 5.94 ± 0.77 km s" 1 for the J = 3 <- 2 
line and uo.54 = —12.00 ±0.32 km s^ 1 and 6V54 = 



8.49±0.81 km s _1 for the J = 5 «- 4 line. Averag- 
ing the results of the two frequency bands, we ob- 
tained the mean line velocity, vo = —12.44 km s _1 , 
and the mean line width, a v = 7.21 km s _1 . We 
then repeated fitting at each pixel of the image 
with the line width and velocity fixed at the mean 
values, a v and vq, solving for the absorption am- 
plitude, AS(v), only. Table Q] lists the measured 
values at the continuum peak pixel (in the map 
center) at each frequency, and the T rot measure- 
ment assuming identical covering factors, f c , for 
the two transitions. If we assume f c = 1, we mea- 
sure T rot = 5.6^0 9 K, and varying f c gives similar 
results. This result is consistent with the predicted 
value of Tom b = 5.14 K and with previo us mea- 
led. 



suremen ts by IWiklind fe Combes 
6.0 K1. ICarilli et al. I (Il 998l T, 
Henkel et al. I (|2009h 



rot 



4.5 



and Muller et al. I (j201ll ) 



T TOt < 

+ 1.5 R N 

Q 6 i±Ji 



However, as shown in Section 01 if the absorption 
does not completely cover the continuum source, 
insufficient spatial resolution can yield biased ro- 
tational temperatures. 

3.2. High Resolution Analysis 

3.2.1. Imaging 

Next, we made images with the higher angu- 
lar resolution afforded by our long interferometer 
baselines. From the line data, we made an im- 
age cube in each frequency channel with 128 x 128 
pixels of 0.05 mas (i.e., for a field of view of 6.4 
mas x 6.4 mas). We adopted identical circular 
restoring beams of 0.65 mas width (FWHM) for 
the 14.5 GHz and 24.1 GHz bands. (The dirty 
beam size was 1.45 mas x 0.46 mas with a posi- 
tion angle of —10 at 14.5 GHz and 1.27 mas x 0.31 
mas with a position angle of —183 at 24.1 GHz.) 

The line data were self-calibrated with the con- 
tinuum at each frequency, and therefore the line 
and continuum maps are accurately aligned with a 
precision that is only limited by the uncertainty in 
the baseline calibration and by the thermal noise 
in the line image (due to the lower signal-to-noise 
ratio for absorption). The bandpass phases are 
calibrated to better than 10 degrees, which trans- 
lates to a positional accuracy of 10/360 of the 
fringe spacing, i.e., ~ 0.03 mas for a ~ 1 mas 
fringe spacing. The positional uncertainty intro- 
duced by the thermal noise in the line image can be 
estimated as half of the beam size #fwhm divided 
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by the signal-to-noise ratio (SNR) of the image. 
With #fwhm = 0.65 mas and SNR^ 6 for absorp- 
tion maps, this results in an accuracy of 0.05 mas. 
For this case, relative positions measured in the 
line and continuum maps are accurate to better 
than 0.06 mas. 

Figure [2] shows the high resolution maps of the 
continuum emission of PKS 1830—211 and the ab- 
sorption in the foreground gas at z = 0.89. The 
continuum maps were obtained by imaging the 
continuum-only data set. The absorption distri- 
bution in each frequency map was obtained by 
plotting absorption depth |A5| measured at each 
pixel of the image (Section I3.2.2|) . Notably, we 
find that the absorption peak position in each fre- 
quency band, in particular at 14.5 GHz, is off- 
set from the continuum peak position. The off- 
sets are significantly larger than the uncertainty 
in the alignment between the line and contin- 
uum maps. In the map of the J = 3 <— 2 
transition observed at 14.5 GHz, the absorption 
peak lies at a positional offset of (Aa cos 6, AS) = 
(0.24 ± 0.04, 0.56 ± 0.05) mas from the continuum 
peak, measured through fitting in the deepest ab- 
sorption channel. For the J = 5 <— 4 transi- 
tion observed at 24.1 GHz, the measured offset 
is (Aa cos 6, AS) = (0.04 ± 0.01, 0.13 ± 0.01) mas 
from the continuum peak to the absorption peak, 
fitted over the averaged map of the 10 channels 
(rs 5 km s _1 width) in which the absorption was 
most prominent. The presence of these offsets is 
a cause for concern that the rotational tempera- 
ture, T rot , measured with lower angular resolution 
may be systematically biased (see Section 0]). 

3.2.2. Spectral fitting 

Since the data of PKS 1830-211 were self- 
calibrated, the absolute positions of both 14.5 GHz 
and 24.1 GHz maps are unknown, and hence we 
need to consider the proper alignment of the two 
frequency maps before we measure T rot . Were 
the continuum emission point-like, or were its 
distribution identical at the two frequencies, self- 
calibration would correctly align the two maps. 
However, the continuum image of PKS 1830—211 
has resolved structure and the relative registration 
of the maps at the different frequencies cannot 
be a priori determined. Based on VLBI images 
of other sources, it is likely that the continuum 
emission of PKS 1830—211 is elongated due to a 



faint one-sided jet emanating from a self-absorbed 
core. Owing to optical depths varying with fre- 
quency for individual components, one expects a 
slight positional offset (along the jet axis) between 
the emission peaks at two frequencies. (Note that 
since absorption depths are small compared to the 
continuum levels, the brightness peak positions of 
the continuum maps are not significantly affected 
by absorption.) 

We measured the absorption amplitude at dif- 
ferent positions across the images, allowing eval- 
uation of different alignments between the two 
images. This was done in the following three 
steps using least-squares fitting of Gaussian pro- 
files to the spectra, in the same manner as for 
the low resolution images. Firstly, we placed the 
measurement position for each frequency map at 
its absorption peak, and solved for the following 
three parameters: the amplitude AS(v) the line 
width (FWHM) a v and the central line velocity 
vo of the absorption. We defined AS(v) as 
AS(v) — Sy me (v) — S c , where Sy lne (v) is the ob- 
served spectrum and S c is the continuum level 
(i.e., AS(v) < for absorption). We then av- 
eraged the results of the two frequency bands to 
determine the mean line width a v and the mean 
velocity vq- We obtained vo — —12.46 km s _1 and 
a v = 6.38 km s^ 1 by averaging i>o,32 = — 12.37 ± 
0.40 km s _1 and a Vl32 = 7.27±0.98 km s" 1 for the 
J = 3 «- 2 line and v 0M = -12.54 ± 0.17 km s" 1 
and (7^,54 = 5.49 ± 0.42 km s -1 for the J = 5 <- 4 
line. Secondly, we repeated the fit at the same 
positions, but fixed the line width and velocity at 
the mean values, a v and Uq> solving for the ab- 
sorption amplitude AS(v) only. Figure [3] shows 
the result of this spectral fit. The obtained values 
are AS 32 (v) = -28.6±2.8 mJy/beam for the J = 
3^2 line and AS 5i (v) = -39.4 ± 2.2 mJy/beam 
for the J = 5 <— 4 line. Finally, using the same 
values of u a bs and a v , we measured the absorption 
amplitude AS(x, y, v) of all individual spectra at 
each pixel position (x, y) across the two images. 
The measured AS(x, y, v) is plotted as absorption 
distributions in Figure [2j 

Previous VLBI observati ons of PKS 1830-211 
at 15 GHz and 43 GHz (iGarrett et al. I Il997 



1998; Jin et al. 112003ft show complicated structure 



of both the SW and the NE images of the ra- 
dio core with extended jet components on a sub- 
milliarcsecond scale. In particular, multi-epoch 
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43 GHz maps by iGarrett et al. I (|l997l) have re- 
vealed that the core has remarkably complex sub- 
structure that changes rapidly with a timescale of 
weeks. The complex morphology of the core sub- 
structure makes it ambiguous to determine the jet 
direction. We therefore determined the jet axis 
based on our continuum maps to lie along the off- 
set between the absorption and continuum peaks, 
and we shifted the position of the T rot measure- 
ment along this axis to allow for the possibility of 
misalignment of the two maps. 

3.2.3. r ro t measurement 

For a given alignment of the two images, T rot 
can be measured at any pixel position, provided 
that both absorption and background contin- 
uum emission have sufficient signal-to-noise ra- 
tios. However, since the detected absorption was 
weak for the J — 3 ■<— 2 line, we fixed the po- 
sition for measuring its optical depth, T32, at its 
absorption peak position in the 14.5 GHz map, 
maximizing the signal-to-noise ratio of the mea- 
surement. We shifted the measurement point of 
T54 along the line shown as a dash-dotted line in 
Figure [2] Both continuum and absorption levels 
were measured at a single pixel, and we performed 
T rot measurements using three different covering 
factors (f c — 1,0.5,0.2), assumed to be identical 
for the two transitions. 

Of course only a very limited range of shifts 
make physical sense, and we required the con- 
tinuum emission to reflect realistic synchrotron 
spectral indices for all pixels with firm detections. 
Our selection criteria required spectral indices, 
— 1 < a < 2 (defined by the power-law spectrum 
S v oc v a ), which yielded four possible alignments 
shown by crosses in the 24.1 GHz map in Figure[2] 
Table [2] lists the measured values from these mea- 
surement points (as well as from excluded mea- 
surement positions with spectral indices outside 
of —1 < a < 2). These alignments yield T rot be- 
tween 1.1 and 2.5 K, which is significantly lower 
than the theoretically-predicted value for the cos- 
mic microwave background, Tcmb = 5.14 K at 
z = 0.88582. (The three excluded measure- 
ment positions in Table [2] near the map center, 
i.e., (|Ax|, |Ay|) < (0.1,0.3) mas, yield spectral 
indices, a, between 2.5 and 2.8, and T rnt values be- 



citation. While this could explain a higher value 
of Tcmb (as in lHenkel et al. II2009I ). it cannot ex- 



tional constant of HC3N, the measured excitation 
temperature may be susceptible to additional con- 
tributions from collisional excitation from (even 
moderate densities of) H 2 gas or local radiative ex- 
explain a 
1. I l2009h . 

plain why our measurement is lower than the ra- 
diation temperature. We discuss the dominant bi- 
ases in the present result in Section |4] By ex- 
ploiting T ro t measurements using Equation [5J we 
assume that the physical properties of the gas is 
the same for the two transitions (i.e., the both 
absorption lines come from the same absorbing 
cloud). We will also discuss this assumption in 
Section HI 

4. BIAS OF TEMPERATURE MEA- 
SUREMENT 

Our high-resolution measurement yields T ro t 
between 1.1 and 2.5 K, which is significantly 
lower than the expected value of Tcmb = 5.14 K 
at z = 0.89. In this section, we discuss several 
possible explanations for this discrepancy. 

4.1. Contribution from non-absorbed con- 
tinuum emission 

Our VLBI maps of PKS 1830-211 have re- 
vealed that the absorption peak positions of the 
HC 3 N(3 <- 2) and HC 3 N(5 <- 4) lines are offset 
from the continuum peak of the SW core. Because 
these two lines were observed contemporaneously, 
temporal variations can be excluded, and these 
offsets indicate that the absorbing cloud of the 
foreground galaxy is covering the core-jet emis- 
sion from PKS 1830—211 only partially, and not 
entirely. 

The synchrotron radiation from the extended 
components of AGNs such as jets is optically 
thin at cm wavelengths, with a spectral index 
a fa —0.7 to —1.0 (defined as S v cx v a ), whereas 
the compact core is generally optically thick and 
synchrotron self-absorption at low fre quencies in- 



verts the spectrum toward a > (e.g., ISchneider 



tween 1.1 and 1.6 K.) We note, as iHenkel et al 



(2009) 3ointed out, that due to the small rota- 



l2006t iRvbicki fc Lightnian1ll979h . Since the core 
is very compact and unresolved on milliarcsecond 
scales, our continuum measurements likely include 
a blend of jet and core emission. 

Blending non-absorbed continuum emission 
with a > within the beam causes the ratio 
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of the optical depths 

T54(^o) _ AS 54 (v Q ) Sc£2 

t 32 (v q ) AS 32 (v ) 5 Ci5 4 

to decrease from the true value, because the con- 
tinuum level 5 Ci 54 relative to S c ,32 (^54 > ^32) is 
increased (due to a > 0), while the absorption 
terms remain unchanged. This biases the value of 
T ro t given by Equation [2] to be lower than the true 
T rot value as we observe. 

In case of our low-resolution measurement (Sec- 
tion [23J 5 we ar e likely more sensitive to extended, 
low surface brightness emission that could have 
a significant contribution of jet emission that is 
also not absorbed, resulting in our measured spec- 
tral index of a = -0.7 ± 0.6 (Table [TJ. The ad- 
dition of non-absorbed continuum emission from 
extended components with a < can falsely in- 
crease the measured T rot . Of course, the rough 
match between the measured and expected T rot 
values would then have to be coincidental. 

4.2. Different absorbing gas 

In our T ro t measurements, we assumed iden- 
tical physical conditions for the gas giving rise 
to the two transitions. Our high resolution maps 
(Figure [2]) show that the absorption peak of the 
J = 5 •<— 4 line is located more closely to the 
continuum peak than the J = 3 <— 2 absorp- 
tion peak (Section 13 . 2 . 1[) . and also the width of 
the absorption line in Figure [3] is narrower for the 
J = 5 4— 4 transition than for the J = 3 <— 2 
transition, which could be interpreted as the two 
absorption lines originating in different gas clouds. 
If this were the case, no relation between the op- 
tical depths of the two absorption lines could be 
assumed, and therefore any value for T rot could 
have resulted from the measurement. Theoreti- 
cally, it is possible to detect only J = 3 <— 2 ab- 
sorption from a given cloud where the J = 5 4— 4 
line is very weak at low temperature, but when 
the J = 5 <— 4 line is detected, the J = 3 <— 2 line 
must also be detected except under highly unusual 
circumstances. Therefore, it is very unlikely that 
we see two different absorbing clouds. 

4.3. Core-shift of PKS 1830 211 and mis- 
alignment at two frequencies 

The observed difference in offset between con- 
tinuum and absorption peaks at the two frequen- 



cies could also be caused by different apparent po- 
sitions of the core at the two frequencies, rather 
than differences in the foreground gas. The radio 
core positions of AGNs are known to shift with fre- 
quency due to frequency-dependent opti cal depths 



of the jets (jBlandford fc Konigl 1 119791 ) . Such a 



core-shift e ffect has bee n obse rved, for example, 
recently by lHada et all (|201lh in M87, but with 



only 0.1 mas shift between frequencies of 15.4 and 
23.8 GHz, comparable to our observing frequen- 
cies. This is well within our resolution, and con- 
sidering the larger distance to PKS 1830—211, this 
effect should be negligible in our data. 

5. CONCLUSIONS 

Our VLBI observations of PKS 1830-211 SW 
have identified offsets between the absorption 
(HC3N J = 3 i- 2 and J = 5 <- 4 lines) and con- 
tinuum peaks (all observed contemporaneously). 
These offsets indicate that absorption in the fore- 
ground galaxy covers the core-jet emission from 
PKS 1830-211 only partially. We have shown 
that insufficient spatial resolution can cause unre- 
solved, non-absorbed continuum emission to con- 
taminate T rot measurements, and that all previ- 
ous T ro t measurements could be affected by this 
bias (although measurements at millimeter wave- 
lengths may be less affected by this bias because 
the jet component becomes less pronounced). Our 
high-resolution measurement yields T rot between 
1.1 and 2.5 K, which is significantly lower than 
the expected value of Tcmb = 5.14 K from the 
cosmic microwave background at z = 0.88582. 
The true rotational temperature of the absorb- 
ing gas is likely higher than this measured T rot 
value, because of blended-in core-jet components 
with spectral index a > 0, which biases the mea- 
surement toward a lower temperature. Our result 
reinforces the importance of high-resolution obser- 
vations of PKS 1830—211 to map the absorbing 
regions in detail, in order to reliably measure the 
rotational temperature and the cosmic microwave 
background. More observations with higher sen- 
sitivity to measure absorption depths as well as 
an absolute-position measurement for correct map 
alignment will further improve the accuracy of the 
measurement. 
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Facilities: VLBA 

APPENDIX 

Our VLBA observations of the southwest (SW) 
core of PKS 1830—211 were used to estimate the 
rotational temperature T rot of the foreground gas 
at z = 0.88582 by measuring absorption depths of 
HC 3 N in the J = 3 <- 2 and J = 5 <- 4 lines. In 
this section, we derive Equation [2] that was used 
to measure the rotational temperature of the gas 
from the observed absorption depths. 

The ratio of the volume densities n of HC3N 
molecules in the lower states of the two transitions 
are given by: 



"-4 
"2 



54 

92 



exp 



-AE 4 
kT Iot 



where gj = 2J+1 are the statistical weights of the 
levels, k is the Boltzmann's constant, and AE 42 
is the energy difference between t he J = 4 and 
2 lev els (AE i2 /k = 3.06 K; iMuller et al~ 
2005). Solving the above equation for the 



J 



2001 



rotational temperature yields: 
-AE 42 



l n (H4 9i 

«2 94 



We now evaluate the density ratio 714/712 in terms 
of the observed optical depths t 32 (V) and t§±{v) of 
the J = 3 <— 2 and J = 5 <— 4 lines, respectively. 
For each transition, the optical depth t{v) is given 
by: 

T W) — \ n[u)dl = k(v)L, 
Jo 

where k{v) is the absorption coefficient, assumed 
uniform, and L is the total path length through 
the absorbing cloud. The absorption coefficient 
k(u) for a transition from the upper level u to the 
lower level I is given by: 



1 cr_ 9v 
,2 



"(/') : -niA ld ( 



8tt v z ul gi 



1 



-hvv.i/kT„ 



where c is the speed of light, v u i is the transition 
frequency, A u i is the Einstein coefficient of spon- 
taneous emission, h is the Planck's constant, T cx 
is the excitation temperature, and 4>(p) is the nor- 
malized line profile, which can be represented for 
any transition as: 



1 



-{u-v y /2a' 



2vrcr' 



where a' v is the frequency half-width at e -0-5 of 
the peak and vq is the line center frequency. This 
can be rewritten in terms of Doppler velocity as: 



1 



-{v-vqY /2a' 



Vq V27TCT 



where a' v is the velocity half-width at e -0 5 of the 
peak, which is related to the full-width at half- 
maximum (FWHM) a v by a v — \/8 In 2 a' v , and vq 
is the line center velocity. Assuming the line pro- 
files of the two transitions are identical in Doppler 
velocity (i.e., a' V 32 = <j' V 5a), then at line center 
(v = vq and v — i/q), the ratio of the line profiles 
is given by the reciprocal ratio of the rest frequen- 
cies: 

032 (vq) _ 

054 (Wo) ^32' 

Therefore, taking the ratio of the optical depths 
T 32(^o) and T54(^o) at line center and assuming 
equal path lengths and excitation temperatures for 
the two transitions yields the following relation: 



7-32(^0) 



V54 \ 93 n 2 g4 A 32 
V32 J 95 n±g 2 A 5i 



1 



-hv 32 /kT., 



\ — Q—hu 5i /hT m 



Solving this equation for the density ratio term 
needed to estimate a rotational temperature, we 
obtain: 



"4.92 



T32(^o) V^32 



93 A 32 
95 A 5i 



1 



-hv 32 /kTe, 



-hu M /UT m 



n 2 gi T)Y2(IA> i \ /':;;> / (/.", A;; \ I 

For J — > J + 1 absorption transitions, the Einstein 
A coefficients are given by: 

64tt 4 



.4 



She 3 



J + l 
2J + 3 



where /1 is the permanent electric dipole moment 
of the molecule, and thus, the ratio of the Einstein 
A coefficients (for J = 2 and J = 4) are given by: 



A 32 
A 54 



3/7 /^32 

5/11 U54 



33/1/32 

35 V^54 



s 



Therefore, we obtain the desired rotational tem- 
perature T r ot in terms of observed absorption opti- 
cal depths of the HC 3 N J = 3^2andJ = 5^4 
lines as: 

Trot = - AE ^ I 

k l n ( 33 §3 T 54 (vo) l- c - h "32/>*T cx \ ' 

\35 g s T 32 (v ) l- e -hv 5i /kT ex j 
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Fig. 1. — Absorption (solid) and continuum (dot- 
ted) maps with a circular beam of 26 mas size 
(FWHM) for the HC 3 N(3 <- 2) line observed at 
14.5 GHz (left) and the HC 3 N(5 <- 4) line at 
24.1 GHz (right). The contours for both contin- 
uum and absorption are spaced by factors of 2, 
starting at 2cr, where a is the r.m.s. noise level 
of the continuum maps, a = 5 mJy/beam. Note 
that the noise level in the continuum map is lim- 
ited by the dynamic range and is comparable to 
the noise level (4 mJy/beam) in the absorption 
map. The peak continuum levels (in the map cen- 
ter) are S c ^2 = 2.066 Jy/beam for the J = 3 2 
line and S c ,$a = 1.429 Jy/beam for the J — 5 <— 4 
line. 




Fast Offsel [mas) Fast Offset (mas) 



Fig. 2. — Same as Figure 1, but with a beam size 
(FWHM) of 0.65 mas, with the same r.m.s. noise 
level of the continuum maps, a — 5 mJy/beam 
(adopted for both absorption and continuum con- 
tours). The peak continuum levels (in the 
map center) are S Ct 32 = 1.340 Jy/beam for the 
,7 = 3^—2 line and iS c ,54 = 0.949 Jy/beam for 
the J = 5 4 line. Measurement positions for 
T Iot are marked by crosses (corresponding to Ta- 
ble [2]), along the dash-dotted line connecting the 
continuum and absorption peaks (see text). 

This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 
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Table 1 

t rot measurement from the low- resolution analysis 



Transition S c 


AS 


a 


fc 


T 


T ro t 


(Jy/beam) 


(Jy/beam) 








(K) 



J = 3 <- 2 2.069 (0.005) 


-0.048 (0.004) 


1.0 
0.5 
0.2 


0.023 (0.002) 
0.047 (0.005) 
0.123 (0.012) 




J = 5 <- 4 1.429 (0.005) 


-0.049 (0.003) 


-0.7 (0.6) 1.0 


0.035 (0.003) 


5 6+ 2 - 5 

°-°-0.9 






0.5 


0.070 (0.005) 


5 6+ 29 

°-°-0.9 






0.2 


0.186 (0.015) 


k q+5.9 

°- y -1.0 



Note. — Measured values from Figured! at the continuum peak pixel in the map center. 
S c is the continuum level, AS* is the absorption level, r is the optical depth of each transition, 
and a is the spectral index derived from the continuum levels at the two frequencies. The 
uncertainties are listed in parentheses. The uncertainty of the spectral index, 0.6, was 
derived assuming a systematic error of 20% in flux scale calibration at each frequency 
band. (Note the flux-scale error does not affect T rot values.) The rotational temperature 
T ro t was measured from the two transitions, assuming identical covering factors f c . 




Velocity (km s L ) Velocity (km s L ) 

Fig. 3. — Baseline-subtracted absorption spec- 
tra (upper panels) and residuals (bottom panels) 
of the HC 3 N(3 <- 2) (left) and HC 3 N(5 <- 4) 
(right) lines toward PKS 1830-211, at the ab- 
sorption peak positions. At this position, the 
continuum levels are S Cl 32 — 0.228 Jy/beam for 
the J = 3 <r— 2 line and S Cy 54 = 0.816 Jy/beam for 
the J = 5 <— 4 line. The fitted spectra are shown 
by smooth curves in the upper panels. 
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Table 2 

t rot measurement from the high- resolution analysis 



Transition 


Aa; 
(mas) 


Ay 

(mas) 


S c 

(Jy/beam) 


AS 
(Jy/beam) 


a 


fc 




T 


Trot 

(K) 


J = 3 <- 2 


0.30 


0.60 


0.229 


-0.029 


— 


1.0 

0.5 


0.13 
0.29 


(0.02) 
(0.05) 
















0.2 


1.01 


(0.26) 


A q+70.6 

%.a _ 12 
6.0t 8 2 2 7 ° 


J = 5 <- 4 


0.20 


0.60 


0.111 


-0.017 


-1.4 


1.0 

0.5 
0.2 


0.17 
0.38 
1.35 


(0.04) 
(0.11) 
(0.53) 




0.15 


0.45 


0.274 


-0.026 


0.4 


1.0 
0.5 
0.2 


0.10 
0.21 
0.67 


(0.02) 
(0.04) 
(0.15) 


o e+0.8 
z -° -0.4 
9 c + 1.2 
z -° -0.4 

9 q + 1.2 

z -° -0.4 




0.10 


0.30 


0.537 


-0.035 


1.7 


1.0 
0.5 
0.2 


0.07 
0.14 
0.40 


(0.01) 
(0.02) 
(0.06) 


1 q+0-3 

i - s -0.2 

1 q+0-3 

-0.2 
1 y+0.4 
1 -' -0.2 




0.05 


0.15 


0.816 


-0.039 


2.5 


1.0 
0.5 
0.2 


0.05 
0.10 
0.28 


(0.01) 
(0.01) 
(0.03) 


1 6 +0 - 2 
i -° -0.1 

-] c+0.2 
1 -° -0.1 

i 4+°- 2 




0.00 


0.00 


0.949 


-0.035 


2.8 


1.0 
0.5 


0.04 
0.08 


(0.01) 

(0.01) 
( c\ ^Q^ 

[U.Uo ) 


1 4+0-2 
x -^-0.1 
1 4+0-2 

i 9+0.2 




-0.05 


-0.15 


0.839 


-0.025 


2.5 


1.0 
0.5 
0.2 


0.03 
0.06 
0.17 


(0.01) 
(0.01) 
(0.03) 


1 3+ - 1 
1 2+ - 1 

ilZ -0.1 
^ -0.1 




-0.10 


-0.30 


0.563 


-0.016 


1.8 


1.0 
0.5 
0.2 


0.03 
0.06 
0.15 


(0.01) 
(0.02) 
(0.04) 


1 9+0.2 
- 1 — ^ -0.1 
i 9+0.2 
ilZ -0.1 
i 1+0.3 
J^- 1 -0.2 




-0.15 


-0.45 


0.289 


-0.011 


0.5 


1.0 
0.5 
0.2 


0.04 
0.08 
0.21 


(0.01) 
(0.03) 
(0.09) 


1 4+0- 3 

-■-•^-0.3 
I 4+0-3 

1 9+0.3 
J^-^ -0.2 




-0.20 


-0.60 


0.115 


-0.009 


-1.3 


1.0 
0.5 
0.2 


0.08 
0.17 
0.52 


(0.04) 
(0.08) 
(0.30) 


2 1+ 12 
2 0+ 17 ' 8 

z - u -0.6 
1 q+60.6 



Note. — Measured values from Figure^ Ax and Ay show the positions (measured in 0.05 mas 
pixel size) of each measurement relative to the continuum peak position in the map center at each 
frequency. S c is the continuum level, AS is the absorption level, r is the optical depth of each 
transition, and a is the spectral index derived from the continuum levels at the two frequencies, 
shown in bold font if it satisfies our selection criteria of — 1 < a < 2. The uncertainties (or noise 
levels) of the continuum S c and the absorption level AS are 0.005 Jy/beam and 0.004 mJy/beam, 
respectively. The uncertainty of each spectral index is 0.6, assuming a systematic error of 20% 
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in flux scale calibration at each frequency band. (Note the flux-scale error does not affect T rot 
values, since T ro t is measured using the relative ratios of the absorption and the continuum levels 
at each band.) We measured the rotational temperature T ro t from the two transitions, assuming 
identical covering factors f c . Measurements for the J = 5 <— 4 transition were made along the 
axis shown in Figure [51 Note that excluded measurement positions with spectral indices outside 
of — 1 < a < 2 are also listed for comparison. 
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